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ABSTRACT 



Context. The Carina Nebula represents one of the most massive star forming regions known in our Galaxy and displays a high level 
of feedback from the large number of very massive stars. While the stellar content is now well known from recent deep X-ray and 
near-infrared surveys, the properties of the clouds remained rather poorly studied until today. 

Aims. By mapping the Carina Nebula complex in the far-infrared, we aim at a comprehensive and detailed characterization of the dust 
and gas clouds in the complex. 

Methods. We used SPIRE and PACS onboard of Herschel to map the full spatial extent (« 5.3 square-degrees) of the clouds in the 
Carina Nebula complex at wavelengths between 70 yum and 500 yum. We use here the 70 fim and 160 //m far-infrared maps to determine 
color temperatures and column densities, and to investigate the global properties of the gas and dust clouds in the complex. 
Results. Our Herschel maps show the far-infrared morphology of the clouds at unprecedented high angular resolution. The clouds 
show a very complex and filamentary structure that is dominated by the radiation and wind feedback from the massive stars. In most 
locations, the column density of the clouds is Nh S 2 x 10 22 cnT 2 (corresponding to visual extinctions of A v S 10 mag); denser cloud 
structures are restricted to the massive cloud west of Tr 14 and the innermost parts of large pillars. Our temperature map shows a clear 
large scale gradient from « 35 - 40 K in the central region to S 20 K at the periphery and in the densest parts of individual pillars. The 
total mass of the clouds seen by Herschel in the central (1 degree radius) region is « 656 000 M Q . We also derive the global spectral 
energy distribution in the mid-infrared to mm wavelength range. A simple radiative transfer model suggests that the total mass of all 
the gas (including a warmer component that is not well traced by Herschel) in the central 1 degree radius region is < 890 000 M s . 
Conclusions. Despite the strong feedback from numerous massive stars and the corresponding cloud dispersal processes that are 
going on since several million years, there are still several 10000M o of cool cloud material present at column-densities sufficient for 
further star formation. Comparison of our total gas mass estimates to molecular cloud masses derived from CO line mapping suggests 
that as much as about 75% of all the gas is in atomic rather than molecular form. 
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1. Introduction 

Most stars form in large clusters or associations, conta ining 
(at least) several thousand stars (e.g., iBriceno et al. 1 I20Q7I) , and 
it has recently become clear that also our sun form ed as part 
of a large star forming complex (e.g., lAdamsll2010l) . In con- 
trast to low-mass star forming regions like Taurus, where the 
individual young stellar objects (YSOs) form more or less in 
isolation and any interaction is minimal, large clusters and 
associations contain high-mass (M > 20 M ) stars. These 
hot and luminous O-ty pe stars profound ly influence their en- 
vironments (see, e.g. Freve r et al.l I2003I) by creating HII re- 
gions, generating wind-blown bubbles, and exploding as super- 
novae. This feedback can disperse the natal molecular clouds 
and thus halt further star formation. However, advancing ion- 
ization fronts and expanding superbubbles can also compress 
nearby clouds and ther eby trigger the fo r mation of new gen- 
erations of stars (e.g., iReach et ail l2004t ICannon et al. 1 12001 



Oev et all [2 005; Preibi sch & Zinneckerlp007t jjpeharveng et a . 



2009; Zavagn o et al.ll2010t iBrand et al 



2011). The interaction 



* The Herschel data described in this paper have been obtained in 
the open time project OTl_tpreibis_l (PI: T. Preibisch). Herschel 
is an ESA space observatory with science instruments provided by 
European-led Principal Investigator consortia and with important par- 
ticipation from NASA. 



of the massive stars with the surrounding clouds and the bal- 
ance between cloud destruction and triggered star formation de- 
termine the charact eristics and the final outcome of the star for- 
mation process (seelDale et al.l2007 t lDale & Bonnelll2008l ; lBatel 
12009b iGritschneder et al.ll2010i for numerical studies). 

At a distan ce of 2.3 kpc, the Cari na Nebula Complex (CNC, 
hereafter; see Smith & Brooks 2008, for a recent review) is the 
nearest star forming region with a large populatio n of very mas - 
sive stars (at least 70 O-type and WR stars; see ISmithl 12006). 
Among these are several of the most massive (M 100 M ) and 
luminous stars known in our Galaxy, e.g., the famous Luminous 
Blue Variable 77 Car, the 02 If* star HD 93129Aa, several 03 
main sequence stars, and four Wolf-Rayet stars. Most of the mas- 
sive stars reside in several loose clusters near the center of the 
complex, including Tr 1 4, 15, and 16, which have ages rangin g 
from <1 to ~ 8 Myr (see lDias et al .120021: IPreibisch et alJ201ld) . 
The CNC contains large amounts of gas and dust clouds, with an 
estimated total mass in the range 350 000 — 670 000 M Q based on 
CO radio maps (e.g. | Grabelskv et al.l 1 9881: lYonekura et al.82 005; 
Smit h & Brooksl2007l) . The very strong radiative and wind feed- 
back from the massive stars has already dispersed much of the 
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original cloud mass in the central region, and drives a large ex- 
panding superbubble, extending over ~ 80 pc (corresponding to 
-2°). 

During the last few years, several surveys have provided a 
wealth of new information on the stella r populations in the CNC. 
The HST survey of the central area (ISmith et al.ll2010aj) and a 
Spitzer survey of the southern parts of the CNC dSmifh et al.l 
2010b) showed that the ionizing radiation from the massive stars 
is currently triggering the formation of a new generation of stars 
in the remaining dense clouds, in particular in the so-called 
"Southern Pillars" region. 

We have performed a very deep near-infrared survey of the 
central ~ 1300 square-arcminute area of the CN C with the near- 
infrared camera HAWK-I at the ESO 8m-VLT dPreibisch et al l 
1201 Id) which is deep enough to reveal all young stars in the 

region with Ay < 35 mag. 

T he Chandra Carina Complex Project (see iTownslev et al.l 
1201 ll for an overview) has mapped the CNC with a mosaic of 
22 individual ACIS-I pointings, using a total observing time of 
1.34 Megaseconds (15.5 days) and covering an area of about 
1.4 square-degrees. With the detection of more than 11 000 X- 
ray emitting young stars, these Chandra data efficiently elimi- 
nate the strong field-star confusion problems that plague visual 
and IR samples and provide, for the first time, a large sam- 
ple of the young stellar population (down to ~ 0.5 M ) in the 
area. A detailed statistical analysis of the X-ray, optical, and 
infrared properties of the detected sources showed that 10714 
of the X-ray source s are very likely young stars in the CNC 



(Bro os et"afll201 lbh . The analysis of the spatial distribution of 
the X-ray detected young stars showed that half of the young 
stellar population resides in one of about 30 clusters or stellar 
group s, while the other half c onstitutes a widely dispersed popu- 
lation (Feigels on et al.l201 ll) . The combination of the X-ray data 
with near- and mid-infrared photometry provided new informa- 
tion about the properties of t he stellar populations i n the entire 
complex dPovich et alJl201 it iPreibisch et al. 201 lbl) and the in- 



dividual clusters Tr 14 and Tr 15 (IWolk et al.ll201 ll IWang et al.l 
1201 lb . These studies showed that the X-ray detected young stel- 
lar populations have ages ranging from < 1 Myr up to « 8 Myr, 
and support the scenario of ongoing, triggered star formation in 
the CNC. 

While these new data sets have strongly boosted the amount 
of available information about the stellar populations, the knowl- 
edge about the clouds in the complex is still much more limited. 
While the IRAS maps and several existing radio maps of CO 
and other molecular lines provide basic information on the large- 
scale morphology of the clouds, only very small areas at the cen- 
ter of the complex have been observed with sub-arcminute an- 
gular r esolution at far-infrared (FIR) and (su b)-mm wavelengths 
so far dBrooks et al.ll2005t iGomez et al1l2010l) . This lack of sen- 
sitive FIR and (sub)-mm data with sufficient spatial resolution 
strongly hampers studies of the cloud structure, and is a very se- 
rious obstacle for investigations of global cloud properties and 
the interaction between massive stars and clouds in the CNC. 

As a first step to improve this situation, we have recently 
used LABOCA at the APEX telescope to obtain a sensitive 
wide-field (1.25° x 1.25°) sub-mm map of the CNC at a wave- 
length of 870//m with 18" angular resolution (= 0.2 pc at the 
distance of the CNC), which provi des the first large-scale survey 
of the dusty clouds in this region (Pre ibisch et al.ll201 lal) . These 
data showed that there are (at least) ~ 60 000 M of dust and gas 
in dense, compact clouds; this represents a large potential for fur- 
ther star formation. While this LABOCA map provided impor- 
tant information about the global properties of the dense clouds 



in the complex, the information that can be retrieved about the 
total cloud phase is limited by two factors: First, the LABOCA 
map is only sensitive to the dense, localized clouds. As a conse- 
quence of the removal of correlated noise in the data reduction, 
structures with angular sizes larger « 2.5' are only partly re- 
covered. Therefore, the more diffuse emission from less dense 
gas is missing from the map. The second limitation is that a 
single-wavelength band map does not allow the determination of 
cloud temperatures. This is critical, because information on the 
cloud temperature is required to compute reliable cloud masses 
from the observed fluxes. The often used approach of simply 
assuming a uniform "typical" cloud temperature (e.g., 15 K) is 
very likely not appropriate in the case of the CNC, where some 
clouds are very strongly irradiated (and thus heated) by numer- 
ous nearby massive stars, while other clouds (especially at the 
periphery of the CNC) experience orders of magnitudes lower 
levels of irradiation (and correspondingly less heating from out- 
side). FIR maps at several different wavelengths can provide 
crucial information about the cloud temperatures and thus al- 
low much more reliable mass and column-density estimates than 
based on single-wavelength data. 

In the present paper we will be discussing observations of 
the CNC performed with the ESA Herschel Space Observatory 
dPilbratt et al.l 120 1 Oh . in particular employing Herschel's large 
telescope and powerful science payload to d o photometry using 
the PACS (Po glitschet al.ll2010h and SPIRE jGriffin et alfc oiO) 
instruments. We will present an overview of the project, deter- 
mine color temperature and column density maps and investigate 
the statistics of these quantities in selected regions of our maps. 
We also investigate global properties of the CNC clouds derived 
from our results. 



2. Herschel far-infrared mapping of the Carina 
Nebula 

In the Herschel open time project OTl_tpreibis_l we used 
PACS and SPIRE to map the entire extent of the CNC in five 
bands between 70 /mi and 500 //m. The observation was per- 
formed on 26 December 2010 in the SPIRE PACS Parallel Mode 
and used 6.9 hours observing time. The J2000 coordinates of the 
aimpoint are RA = 10h44m03s, Dec = -59d30m00s. The fast 
scan speed (60 arcsec per second) was used to map a « 2.3°x2.3° 
area[3, which corresponds to a physical region of « 92 pc x 92 pc 
at the distance of the CNC and covers the full spatial extent of the 
CNC. The mapped area als o includes the HII region Gu m 31 (= 
Ced 108 = NGC 2599; see ICederbladll 1 9461: iGumll 1 955l) around 
the young stellar cluster NGC 3324 at the north-western edge of 
the Carina Nebula. 

The observation was done with two scan maps of the same 
area, one with nominal scan direction, and the other with or- 
thogonal scan direction, in order to remove more efficiently the 
stripping effects due to the 1 // noise and to get better coverage 
redundancy. 

We used the 70 //m PACS channel as the blue band, in or- 
der to get the widest possible wavelength range. The wave- 
length bands of the five maps resulting from our observations are 
[60 - 85] /mi, [ 1 30 - 2 10] //m, [210 - 290] /mi, [300 - 400] //m, 



1 Note that due to the fixed 21' separation of the SPIRE and PACS 
focal plane footprints on the sky, the areas observed by each individual 
instrument are somewhat larger and show some offset with respect to 
each other; the mentioned w 2.3° x 2.3° area is the region that is covered 
by both SPIRE and PACS . 
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and [420 - 610] /mi, and will be denoted as the 70 //m, 160 /mi, 
250 //m, 350 //m, and 500 /mi band in the following text. 

T he data reduction was performed with the HIPE v7.0 dOttl 
1201 Oh and SCANAMORPHOS vlO.O (iRoussell I20I2I) software 
package^. From level 0.5 to 1 the PACS data was reduced using 
the Ll_scanMapMadMap script in the photometry pipeline in 
HIPE with the version 32 calibration tree. The level 2 maps were 
produced with SCANAMORPHOS with standard options for par- 
allel mode observations and the galactic option to preserve 
brightness gradients over the field. The pixel sizes of the two 
PACS maps were chosen as 3.2" (for the 70 /mi map) and 4 .5" 
(for the 160 /mi map), as suggested in Trafic ante et al. d201 ll) . 

The level SPIRE data were reduced with an adapted ver- 
sion of the HIPE script rosette_obsidl&2_scriptJevell included 
in the SCANAMORPHOS package. The final maps were pro- 
duced by SCANAMORPHOS with standard options for parallel 
mode observations and the galactic option as for the PACS 
data. The pixel sizes were chosen as 6", 8" and 11.5", for the 
250, 350, and 500 //m band, respectively. 

The quality of the final maps turned out to be very good. 
In order to characterize the point-spread-function (PSF) in the 
maps, we determined the full-width at half-maximum (FWHM) 
for a number of isolated point-like sources and found values of 
« 10" for the 70 /mi map, a 15" for the 160 /im map, « 20" 
for the 250 /mi map, * 26" for the 350 //m map, and » 36" for 
the 500 //m map. These values show that the image quality of the 
PACS maps is slightly worse than the telescope diffraction limit 
(as expected for observations performed in fast scan mode), but 
the SPIRE maps are close to the ideal quality. With their un- 
precedented angular resolution of ~ 10" - 36", corresponding to 
physical scales of ~ 0.1 - 0.4 pc, these Herschel maps represent 
the first detailed, deep, and spatially complete FIR maps of the 
fullCNC. 

The levels of background cloud fluctuations, determined as 
the standard deviation of the pixel value in apparently empty 
regions in the maps, are found to be w 2.3 mJy/px at 70/jm, 
«: 5.0 mJy/px at 160 /mi, w 84 mJy/beam at 250/mi, « 
83 mJy/beam at 350 /mi, and « 73 mJy/beam at 500 /mi. 

A potential problem in the calibration of the Herschel maps 
is the possible presence of zero-level offsets (see Bernard et al.l 
l2OT0h . No correction for these offsets is currently implemented 
in the Herschel data reduction tools. There are several reasons 
why we can expect the uncertainties related to offsets to be only 
moderate for our specific data set and the results obtained in this 
paper. One reason is that by the use of the the galactic op- 
tion in SCANAMORPHOS the baseline subtraction is modified 
with the aim to preserve large-scale structures as well as possi- 
ble. Another important aspect is that our maps cover the entire 
spatial extent of the clouds associated with the Carina Nebula. 
The periphery of our maps shows low-intensity galactic back- 
ground emission that is not related to the Carina Nebula and 
thus not relevant for the aims of our study. Therefore, the emis- 
sion from the relevant clouds associated to the Carina Nebula 
can be well separated from the galactic background in our maps. 
We also compared our Herschel 70 /mi map to the IRAS 60 /mi 
map and found very good agreement (within about 1%) of the 
fluxes integrated over large regions. Finally, we note that in our 
temperature and column density determination described in the 
next section even a 20% error in the intensity causes just a ~ 5% 



2 see http : //www2 . iap . fr/users/roussel/herschel/ 

3 Since the Carina Nebula is close to the galactic plane (b 0.6°) 

and near the tangent point of a spiral arm, the background consists of 
numerous distant cloud complexes. 



change in the derived color temperatures and a ~ 12% change in 
the derived cloud column densities. 

As soon as the Planck data get available, a more detailed 
analysis of the possible effects of offset corrections will be per- 
formed. 



3. Cloud morphology 

The complete set of all five Herschel maps resulting from our 
observation is shown in Figures Al to A5 in the appendix. In 
Figure [1] we compare the 70 pm and 500 pm Herschel maps to 
an optical image and the Spitzer 5.8 /mi image. The remarkable 
similarity between the Herschel 70 /mi image and the Spitzer 
5.8 pm image demonstrates that most of the clouds are only mod- 
erately dense, i.e., are already transparent in the Spitzer 5.8 /mi 
band. This point will be discussed in more detail in Sect. 5. 

In Fig. [2] we show an RGB color composite of the optical, 
70 /mi and 160 /mi emission. This illustrates the spatial rela- 
tion between the hot (~ 10 4 K) Ho- emitting gas and the dense 
cool/cold gas. It shows how the hot gas fills the interior of the 
bubbles in the cloud structure. 

The comparison of the Hersche l 70 um maps to our 
LABOCA 870 /an map presented in IPreibisch etail (1201 lab 
shows that the detectable sub-mm emission is restricted to a 
small volume fraction of the cloud shown by Herschel. The vol- 
ume filling factor of the dense clouds is very low, and this dense 
material is highly fragmented and dispersed throughout the area. 
Herschel reveals the much more widespread diffuse gas at lower 
densities. 

Figure [3] is an unsharp-mask filtered version of the 70 /mi 
map that highlights the small-scale structure of the clouds. It 
clearly shows that the clouds are highly filamentary and have a 
very complex morphology. Herschel mapping of other star form- 
ing regions has shown that molecula r clouds generally exhibit an 
extensive filamentary structure (e.g. lAndre et"aDl2010l) and that 
the typical width of these filaments is of the order of ~ 0.1 pc 
(lArzoumanian et aU201 ll) . This implies that the width of the in- 
dividual filaments can be only marginally resolved in our 70 /mi 
map and remains unresolved in the longer wavelength maps of 
the CNC. 

Based on the observed global cloud morphology, we define 
several sub-regions of the CNC and discuss them separately. 



3.1. Central region 

In optical images of the Carina Nebula, the "V"-shaped dark 
clouds just below 77 Car constitute a very prominent feature. The 
Herschel image shows that these dark clouds are in fact of only 
moderate density. As we will find below, their typical visual ex- 
tinction is just Ay ~ 3 - 5 mag. 

The nebulosity to the north of r\ Car is the well known 
"Keyhole Nebula", which was first described by J. Herschel in 
1840 (a reproduction of Herschel's drawing an d a comparison to 
a mod ern photograph can be found as Fig. 4 in Smit h & Brooksl 
2008). While this nebula contains large amounts of hot 
gas, producing c opious Hg emission, it also contains cold 
molecular clouds JCox & Bronfmanlll995t iBrooks et all 12005; 
iGomez et al1l2010t IPreibisch et al.ll201 lal) 

To the south of 77 Car, the Herschel image shows a region of 
low emission. This seems to represent a bubble of lower density, 
that was probably created by the wind- and radiation-feedback 
from the numerous very massive stars in the Tr 16 cluster. 
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Fig. 1. Optical, near-, and far-infrared view of the Carina Nebula. All four images show the same field of view of 
3° x 3°. North is up and east to the left. Upper left: Optical image obtained from the ESO Photo Release 1145 
(http://www.eso.org/public/images/eso0905b/; ESO/Digitized Sky Survey 2, Davide De Martin). Upper right: 5.8 pm image con- 
structed by us from Spitzer data obtained from the public archive. Lower left: Herschel 70 pm map. Lower right: Herschel 500 yum 
map. The Spitzer and Herschel images are displayed with a square-root intensity scale. 



3.2. Southern Pillars (SP) 

The cloud structure in the so -called "Southern Pi llars" (SP) was 
already discussed in detail bv lSmith etall (1201 Obi) on the basis of 
the Spitzer maps. The cloud morphology we see in our Herschel 
maps is very similar. 



it as the "Northern Cloud" (NC). The eastern edge, where this 
cloud faces the stellar clusters, is the site of a prominent Photon 
Domi nated Region (PD R), which has been studie d in some detail 
bv iBrooks etaU (120031) and lKramer et al.1 J2008). Our Herschel 
maps show that the cloud extends (at least) about 15' (~ 10 pc) 
to the west. 



3.3. Northern Cloud (NC) 

The cloud to the west of the center, close to the stellar cluster 
Tr 14, is clearly the den sest and most massive cloud structure in 
the complex. Following ISchneider & Brooksl (120041) . we denote 



3.4. Bubbles north ofi] Car (NB) 

The field north of rj Car also shows a few prominent pillars, but 
the general structure is dominated by a web of bubbles and arcs. 
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Fig. 2. RGB composite of the red optical (DSS) image in blue, the Herschel 70 /mi image in green, and the Herschel 160 fim image 
in red. The image shows a 2.9° x 2.6° field of vie w. North is up and east to the left. The green plus signs mark the positions of the 
high-mass stellar members of the CNC as listed in ISmithl(l2006h . 



We will therefore denote this region as the "Northern Bubbles" 
(NB). The diffuse appearance of the clouds and the lack of sub- 
mm emission from dense clumps in this region suggests that 
the density contrast of the clouds is smaller than in the SP re- 
gion, where numerous dense and massive pillars are present. 
This morphological differences may indicate that the massive 
star feedback in this area is different from the strong irradiation 
that shapes the surfaces of the southern pillars. 

The bubble-like structure of the clouds is partly related to in- 
dividual HII regions that have swept up shells of dense clouds 
around them. The two most prominent HII regions in this area 
are the extended nebula surrounding the Wolf-Rayet star WR 23 
(=HD 92809; spectral type WC6), and the RCW 52 (=Gum 32) 
nebula (which is illuminated by the 07V star LSS 1887). Both 
regions are clearly surrounded by bubbles of gas and dust, which 
have been interpreted as being mainly created by the action 



of the stellar winds of t he massive stars on their environment 
( Cano a et al.ll2005ll201 lh . 



The irregular filamentary structure could be caused by 
evolved bubbles that already broke up into pieces. We note that 
this region contains the c luster Tr 15 and Bo 10, which have esti- 
mated ages of 7 - 8 Myr dDias et al.ll2002b iPreibisch et al.l201 id) 
and are thus several Myr older than the clusters Tr 14 and Tr 16 
in the center of the Carina Nebula. This could suggest that at this 
older age, the stellar-wind feedback (especially from the evolved 
massive stars) may play an important role. However, the validity 
of this interpretation remains unclear, since recent studies sug- 
gest that massive star feedback is usu ally dominated by ioniz- 
ing radiation (e.g. Marti ns et al.ll20l6T) and the effects of ste llar 
winds are of secondary importance (e.g. Martins et al. 2012|). 
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Fig. 3. Unsharp-masked version of the Herschel 70fim image with a logarithmic intensity scale image. The different regions men- 
tioned in Section 3 are marked by the black boxes. The positions of rj Car, the clusters Tr 14, Tr 14, and Bo 10, as well as the HII 
regions around WR 23 and RCW 52 are marked. North is up and east to the left. 



3.5. The large southern bubble 

In the Spitzer images, the large elongated bubble south-west of 
rj Car is a quite prominent hole structure. Our Herschel maps 
confirm that this structure is truly an almost empty bubble, and 
not caused by the shadow of a dark cloud. Be sides the small 
dense i nfrared dark cloud described already in Preibisch et al.l 
( 1201 lab , only very low levels of emission are seen in the in- 
ner region of this bubble in the Herschel maps. The Herschel 
maps suggest that the density increases again at the southern 
edge, i.e. the bubble wall may be closed at the southern edge. 
The diameter along the major axis is about 50', corresponding to 
a physical length of about 33 pc. 

The comparison of our FIR maps to the optical image (Fig. 2) 
shows that the Ha emission is well aligned with the inner edges 
of the inner cloud surface; as the optical emission is strongest at 
the edges and lower in the central area, this suggests a shell-like 
morphology for the hot gas. 

The origin of this bubble remains unclear. Our search of the 
SIMBAD database revealed only three known early-type stars 
inside this bubble: the 08 star HD 305438, the Bl/2 II star 
HD 92741, and the B2 III star HD 92877. None of these seems 
to be massive enough to have created this large bubble. One pos- 



sible explanation is that the hot gas from the winds of the nu- 
merous massive stars in the central area is leaking out and filling 
this bubble. An alternative possibility is that this bubble repre- 
sents the southern lobe of a large bipolar HII region that is driven 
by the massive stars in Tr 16. 

3.6. The Gum 31 bubble 

In the north-western part of our Herschel maps, the prominent 
circular bubble around the HII region Gum 31, which is pro- 
duced by the young stellar cluster NGC 3324, is clearly the 
dominating cloud structure. It represents a nice example of a 
"perfect bubble" of den se clouds around an HII region (see, e.g., 
Deha rveng et al.ll2009l for other examples) that has been swept 
up by the expansion of the HII region. 

It is interesting to note that the Herschel maps show elon- 
gated filamentary clouds that seem to connect this bubble to the 
central regions of the Carina Nebula (see Fig. 2). The idea that 
the Gum 3 1 bubble is part of the giant Carina Nebula cloud com- 
plex is supported by the available information about the radial 
ve locities of the molecu lar clouds. According to the CO maps 
of lYonekura et al.l (|2005), the clouds surrounding Gum 31 have 
radial velocities in the range Vlsr = -24 ... - 21 km s _1 , which 
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is very similar to that of the clouds in the central regions of the 
Carina Nebula, for which a range of Vlsr = -26 ... - 18 km s 
is found. Another interesting aspect is that the densest and most 
massive parts of this bubble are found at the south and south- 
eastern regions of the rim, i.e., just along a line connecting the 
center of the bubble to the center of the Carina Nebula. This may 
indicate some kind of interaction between the clouds associated 
with the Carina Nebula and the Gum 31 bubble. A more com- 
prehensive and detailed investigation of these aspects will be the 
topic of a forthcoming paper. 



4. Cloud colors and temperatures 

4.1. Morphology in multi-color FIR images 

Figure |4] shows an RGB composite of the central region of the 
CNC created from our three shortest Herschel wavelengths, in 
which the colors trace the different cloud temperaturesj. The 
bluish emission traces the relatively warm gas in the Keyhole 
Nebula Region, and also in front of the eastern edge of the 
Northern Cloud, which is strongly irradiated by the stars in the 
cluster Tr 14. 

The red tones in this image trace the coldest cloud mate- 
rial, which is concentrated in the densest parts of the base of 
the pillars in the Southern Pillar region, and also in a number of 
compact clouds north of the Northern cloud. 

There are several reddish dense cold clouds to the north-west 
of center; in contrast to similarly dense clouds in the Southern 
Pillar region, these do not show a pillar-like structure. 

4.2. Color temperatures 

In order to estimate the local temperature of the clouds in our 
map, we determined a color temperature from the ratio of the 
observed fluxes in two Herschel bands. Due to the high level of 
radiative feedback from massive stars, which heats the clouds, 
and the fact that most clouds have only moderate column den- 
sities, we can expect that most of the clouds in the CNC should 
be not extremely cold (such as infrared dark clouds which of- 
ten show temperatures of T < 15 K), but somewhat warmer. We 
will therefore employ the ratio of the 70 fim versus 160 //m map 
intensities, which is well suited to measure temperatures in the 
range ~ 20-75 K0- This choice has the advantage that the result- 
ing color temperature maps provide the best angular resolution 
that can be obtained from our Herschel maps. 

In order to construct a flux-ratio map, we first convolved the 
70 fim map with a kernel as discussed in lAniano et al. I (l20Tlt) 
to match the angular resolution of the 160 fim map. Following 
the procedure for the temperature dete rminations in the Hi-GAL 
data described bv lBernard et al.l(l2010l) . we divided the fluxes de- 
rived from the maps by the recommended color-correction fac- 
tors of 1 .05 and 1 .29 for the 70 jum and 160 fim PACS bands that 
were derived from a comparison of PACS data to Planck data. 



4 Evaluating the peak of the Planck function B V (T) shows that the 
70yum band is most sensitive to cloud temperatures of T ~ 72 K, the 
160 fim band to T ~ 32 K, the 250yum band to T ~ 20 K, the 350^111 
band to T ~ 15 K, and the 500 yum band to T ~ 10 K. 

5 For a dust dust emissivity index of ft = -2.0, the 70/um versus 
160 fim flux ratio is 1 for T = 28.3 K, 10 for T = 74.8 K, and 0.1 for 
T = 18.0 K. 



Assuming that we see optically thitf] thermal dust emission, 
the ratio of the observed fluxes can be related to the temperature 
via the formula 

WOpm) K B Y {T){10fim)K v {10fim) = 
7 y (160jum) 00 B v (r)(160//m)xv(160/an) 1 

We further assume that the dust emissivity follows a power- 
law (k v oc The value of the dust emissivity index /3 is not 
very well known and under debate; typical observationally de- 
termined values range from ~ 1 to ~ 3, an d it is possible tha t 
P is anti-correlated to the temperature (e.g., IShettv etldl 12009). 
Here we assume f3 = 2, what should be appropriate for the rela- 
tively warm clouds in the Carina Nebula and has been confirmed 
for the not too cold (T > 20 K) clouds in HII regions (see, e.g., 
lAnderson et al.ll2Q10h . Equation Q] allows to determine the tem- 
perature in each pixel by comparing the observed flux ratio with 
a pre-computed table of flux-ratio values on a finely space tem- 
perature grid. 

A complication of this temperature determination method re- 
sults from the noise in the Herschel maps. In regions with low 
surface brightness, statistical fluctuations can lead to large errors 
of the derived temperatures due to division by very small num- 
bers. Furthermore, some outer regions of the maps contain pix- 
els with negative values, preventing any temperature determina- 
tion. In order to avoid these problems, we restricted the compu- 
tation of color temperatures to pixels with intensities above the 
limits of 7 v (70yum) > 0.001 Jy/square-arcsec and 7 y (160jum) > 
0.01 Jy/square-arcsec. These limits are several times above the 
noise-level in the maps and also serve to separate the emission 
of the clouds in the CNC from the large-scale galactic back- 
ground emission. This minimum intensity condition is fulfilled 
for 3 000757 pixels in our map, i.e. we can derive color temper- 
atures for a total area of about 2.37 square-degrees (« 44% of 
the full extent of our maps). Since we are only interested here 
in the properties of the relatively bright clouds associated with 
the Carina Nebula, these temperature maps cover nevertheless 
almost the entire area of interest. In the much fainter regions, 
predominantly around the periphery of our map, for which we 
cannot determine temperatures, the emission is dominated by 
the galactic background, which is not of interest for our current 
study. 

A potential problem for the color temperatures may be ex- 
cess emission from very small dust grains. The absorption of 
individual optical or UV photons from the surrounding stellar 
radiation field can cause strong temperature excursions of the 
smallest grains (with sizes below about 10 nm). Due to these 
stochastic temperature fluctuations by single-photon heating, the 
very small grains may not reach an equilibrium temperature. 
The emission from very small grains in the high-temperature 
state, immediately after a photon absorption, can lead to an ex- 
cess of mid-infrared radiation in the total emission spectrum of a 
cloud with a distrib ution of grain sizes. As discussed in detail in 
iDraine &"Q (120071) . two mechanisms play a role. First, photon- 
heated very small grains consisting of polycyclic aromatic hy- 
drocarbon (PAH) can emit strong mid-infrared band radiation, 
in particular around the wavelengths of 11.3, 12.7, and 17 fim. 
Second, very small carbon and silicate grains may cause an ex- 
cess of mid-infrared continuum emission. Whereas the Herschel 
bands do not contain strong PAH emission features, they may be 
affected by continuum excesses. The strength of excesses from 
very small grains depends on the ambient radiation field. When 

6 The assumption that the clouds are optically thin will be justified 
and confirmed in our analysis of the derived column densities in Sect. [5] 
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Fig. 4. RGB composite of the central region (1.40° x 1.43°; north is up and east to the left) constructed from our Herschel 70/im 
(blue), 160 fim (green) and 350 /urn (red). 



the rate of starlight heating is large (as is the case in the Carina 
Nebula), the relative amplitude of the temperature fluctuations of 
the grains decreases and the the steady state equilibrium temper- 
ature approximation can beco me valid even for ve ry small grains 
(see discu ssion in Sect . 3 of [Draine "&LHl2007t) . According to 
Table 4 in Draine & Li (2007) and in the presence of a radiation 
field that is 100 times stronger than in the local solar environ- 
ment, a variation of the abundance of very small grains by a 
factor of 10 leads to variations in the emissivities of 17% for the 
70yum PACS band, 7% for the 160 /mi PACS band, and < 5% 
for the SPIRE bands. We note that a 20% error in the 70 /mi in- 
tensity causes a ~ 5% change in the derived color temperature. 
This represents a lower limit to the accuracy of the temperature 
determinations. 

Figure [5]shows a map of the resulting color temperature val- 
ues. A histogram of the derived color temperatures can be found 
in Fig. [7] The color temperatures range from 14 K to 83 K; the 
mode and the median of the distribution is 29 K and 27.3 K, 
respectively. Only 4% of the pixels have T < 20 K, and the 
fraction of pixels with temperatures of T > 40 K [50 K] is just 
0.39% [0.01%]. 

In order to estimate the uncertainties of the derived color 
temperatures, we first note that the flux calibration of the 
Herschel maps is thought to be accurate to ~ 10 - 20%; at 
T ~ 30 K, a 20% intensity error causes a 5% (~ 1.4 K) change of 



the derived color temperature. Besides the above mentioned un- 
certainties related to possible excess emission from very small 
grains, another contribution results from the photometric cali- 
bration uncertainties of the Herschel bands. For temperatures 
below ~ 20 K, the quickly rising photometric color correction 
terms cause higher uncertainties and a possible bias towards 
over-estimated temperatures; however, since only 4% of the pix- 
els in our maps yield such low temperatures, this should not be 
a serious problem in the analysis of the global properties of the 
clouds presented in this paper. Furthermore, since our color tem- 
peratures are derived from the Herschel maps with the shortest 
wavelengths, they may be biased towards higher dust temper- 
atures (and consequently lower column densities). Preliminary 
results from our more detailed investigation (that will be pre- 
sented in a subsequent publication) suggest that any such bias 
is rather small. In summary, we estimate the uncertainties of the 
derived temperatures to be about ~ 10 - 15%. 



We also would like to point out that the temperature values 
we determined represent averages over the line-of-sight through 
the full depth of the clouds. Since in some locations there may be 
several individual clouds that are seen projected onto each other, 
the physical meaning of the derived temperature is limited, but 
this is a general problem for the study of cloud properties, and 
not specific to our data set. 
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Fig. 5. Map of 70-160 /im color- temperatures. Pixels with insufficient 70 /-im or 160 jum intensities for a reliable determination of 
the intensity ratio are left white. The blue plus signs mark the locations of the high-mass stellar members of the CNC as listed in 
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Fig. 6. Map of column densities (A^h) with symbols as in Fig. [5] Pixels with insufficient 70 /im or 160 fj.m intensities for a reliable 
determination of the intensity ratio are left white. 



9 



T. Preibisch et al.: Herschel far-infrared observations of the Carina Nebula complex 



Our temperature map in (Fig. |5) shows a clear systematic 
temperature gradient from the central regions (where most of the 
hot, massive stars are located) to the periphery of the complex. 
A similar, but much steeper temperature gradient is seen within 
the boundaries of the Gum 3 1 nebula. This reflects the fact that 
the number of high-mass stars in NGC 3324 (a bout 3) is much 
smaller than the 126 O-, early B-, and WR stars (ISmifhll2006l) in 
the central regions of the CNC. 

5. Cloud column densities 

With the temperature estimate we can now proceed and com- 
pute the column densities of the clouds. The observed surface 
brightness in our maps can be converted to the beam-averaged 
hydrogen column density via the formula 

^Vh — 2 A^h 2 — 2 F " R , (2) 

By(T d ) LlK v nm n 

where Q is the beam solid angle, R is the gas-to-dust mass ra- 
tion (assumed here to be 100), K v is the dust opacity, and /i is the 
mean molecular weight. For consiste ncy with our pre vious anal- 
ysis of the LABOCA sub-mm data (Preibis ch et al.ll201 Tab , we 
use the Ossenkopf & Henning] ( 1 19941) dust model for grains with 
thin ice-mantles that coagulated at a density of n — 10 6 cm~ 3 
which gives a dust opacity of /e v (160 /mi) = 40.5 cm 2 /g. We note 
that, due to the unknown details of the chemical composition 
and physical structure of the dust grains, opacity values suffer 
from uncertainties of (at least) a factor of ~ 2, resulting from the 
dependence of the dust opacity on the details of the grain proper- 
ties. This causes uncertainties by a factor of ~ 2 for the derived 
cloud masses. 

Figure [6] shows the resulting map of cloud column densi- 
ties. We assume that column density and visual extinction are 
related via the "canonical" relation Ay = 1 mag <-> Nh 
■ X 1 21 cm' 2 <-> N H , = 1 x 10 21 cirT 2 (see, e.g.. lBohlin et alj 
Il978l) . For the analysis of the statistics of the column densi- 
ties, we excluded a 1' X 1' region centered on r\ Car, which is 
a source of strong FIR emission. Fig. Qshows the histogram of 
the column density values. We find the mode of the distribution 
at Ay =1.0 mag and a median value of Ay = 1.17 mag. The 
distribution shows a second peak near Ay « 2 mag, and drops 
steeply towards higher column densities. In the A v ~ 3 - 30 mag 
range, this drops can be well described by a power-law relation 
of the form 



Table 1. Integrated cloud masses for different extinction thresh- 
olds. 



oc (logCJVn))" 



(3) 



d log(AQ 

d l0g(/VH) 

with an exponent a w -2.9. This slope is conside r ably s teeper 
than column density slopes derived by iHill et all (1201 ll) from 
He rschel observation s of the Vela C cloud complex (a w -2.0) 
or iKainulaine n et alj d201 ll) from an extinction map of the 
Ophiuchus cloud (a « —1.5), Although a direct comparison of 
different star forming regions is difficult due to the unavoidable 
differences in the observations and data analysis, this result sug- 
gests that the CNC contains a smaller fraction of the total cloud 
mass at high densities than the other mentioned cloud com- 
plexes. This may be the consequence of the fact that the CNC is 
already a somewhat evolved star forming region (star formation 
has started already about 8 Myr ago), much of the original dense 
cloud mass has already been dispersed by the very high level 
of massive star feedback, and the today present dense clouds are 
less the remnants of the original cloud but more the result of stel- 
lar feedback concentrating and compressing some fraction of the 
remaining clouds. 
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Mm 
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32 300 
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NB 
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4075 
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186 700 


46400 


7325 



Our finding that just 1% of the pixels have Ay > 5.26 mag, 
and only 0.1% have Ay > 12 mag agrees with our estimate of 
the (generally moderate) cloud extinction mentioned above and 
can also be used to provide an (a-posteriori) justification of the 
assumption of optically thin FIR emissi on. The 70 urn opacity of 
the grains in the above mentioned Ossenkop f & Henningl (ri994) 
dust model shows that the optical depth of 99.9% of the pixels 
in our maps (which have Ay < 12 mag) is r(70yum) < 0.036. 



6. Column-density- and temperature-distributions 
and integrated cloud masses in the different 
parts of the complex 

We can now compare the properties of the cloud structure in 
the different parts of the CNC as defined in Sect. 3. First, we 
determine the total cloud masses in each region by integrating 
our column densities over the corresponding area. For the total 
cloud complex associated to the Carina Nebula we define a cir- 
cular region with a radius of 1 degree, centered on the position 
R.A.=10h44ml6s, DEC=-59d38m28s (this point is 6.5' north- 
east of rj Car). This region includes nearly all the clouds in our 
map, but excludes the clouds around Gum 3 1 . Integrating the de- 
rived column densities over this circle yields a total cloud mass 
of M tot = 655 700 M Q . For the Southern Pillars region we find 
209 900 M , and for the Northern Cloud 92 8OOM . The inte- 
grated mass in the Northern Bubbles region is 161 500 M Q , and 
for the Gum 31 region we derive 186 700 M . 

While the derived integrated masses should be largely free 
from background contamination (due to the flux thresholds we 
applied in our analysis), the gas shows a wide range of densities, 
from very diffuse and extended low-density filaments to quite 
massive clouds (like the head of the Northern Cloud). Numerous 
recent investigations of the relation between cloud structure and 
the star formation process lead to suggestions of a (column) - 
density threshold c louds have to exc eed in order to allow ac- 
tive star formation. Lad a et al.l (1201 Oh found that the star forma- 
tion rate in molecular clouds is linearly proportional to the cloud 
mass above an extinction threshold of Ay w 7 mag (correspond- 
ing to a gas surface density thres hold of ^ 100M o pc~ 2 ). The 
study of Kainulainen et alj d201 ll) suggested the transition be- 
tween diffuse clouds and bound clouds (which may be the site 
of stars formation in the near future) near Ay * 3 mag. We have 
therefore determined the integrated masses in the above defined 
regions for all pixels exceeding the column density thresholds of 
Ay = 3 mag and Ay = 7 mag. These values are listed in Tab.Q] 

For the full CNC region, the fraction of mass above Ay = 
3 [7] mag is 16% [3.4%]. While these fractions are very sim- 
ilar for the Southern Pillars and just slightly higher for the 
Gum 31 clouds, they are considerably higher (48% [15%]) for 
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Fig. 7. Histogram of 70-160 fim color temperature values (left) and column density values (right) derived from our maps. 
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Fig. 8. Histogram of 70-160 fim color temperature (left) column density values (right) for the individual subregions defined in 
Sect. 3. 



the Northern Cloud and considerably lower for the Northern 
Bubbles (2.5% [0.1%]). 

In Fig.[8]we compare the distributions of pixel temperatures 
and column densities for the four different regions. This shows 
considerable and remarkable differences between these regions. 
The clouds in the Southern Pillars show relatively high cloud 
temperatures. We interprete this as a signature of the strong ir- 
radiation that heats these clouds. In combination with the rela- 
tively low column densities this suggests that a large fraction of 
this cloud material is currently in the process of being evaporated 
due to this heating. 



The Northern Bubbles are characterized by systematically 
lower cloud temperatures. The distribution of column densities 
is quite similar to that in the Southern Pillars for low extinction 
values (Ay 5 2 mag), but drops off very quickly for higher val- 
ues. This implies that the density contrast in this region is con- 
siderably smaller than in the Southern Pillars. This is probably a 
consequence of the quite different levels and nature of massive 
star feedback imposed onto these clouds. In the Southern Pillars, 
the clouds are directly irradiated, what heats their surfaces and 
can lead to considerable cloud compression in some locations. 
The Northern Bubbles, on the other hand, are apparently in a sit- 
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uation of less strong irradiation and their structure seems to be 
dominated by matter flows. 

The temperature distribution in the Northern Cloud is quite 
similar to that in the Southern Pillars. The same is true for the 
distribution of column densities for values up to log (Nu) < 
22.3. For higher densities, however, the Northern Cloud shows 
a prominent excess peak near log(A^e) < 22.4, which can be 
associated to the dense massive center of this cloud. 

Finally, the clouds around Gum 3 1 show relatively low tem- 
peratures (with a mode value of 21 K), and intermediate column 
densities. 

It is also interesting to compare the masses derived here 
from our Herschel maps to those determined from our previ- 
ous LABOCA map. The total mass of all clouds for which our 
LABOCA map of the c entral 1.25° X 1.25° o bserved region was 
sensitive is 60000M o dPreibisch et al.ll201 lah . Considering the 
same area in our Herschel maps, we find here a total mass of 
427 000 Mq, with 90200M Q above the A v > 3 mag threshold. 
This agrees well with the expectation that LABOCA traces only 
the mass of the denser clouds. 



7. Total cloud mass derived from a simple radiative 
transfer modeling 

If we want to determine the total mass of all the clouds in the 
CNC, we have to take into account that Herschel is not very sen- 
sitive to diffuse warm gas at temperatures ^ 80 K. However, due 
to the very high level of massive star feedback, a substantial frac- 
tion of the total gas mass could be at such relatively high temper- 
atures. This gas should then be in atomic rather than molecular 
form. In order to estimate the total mass of the clouds, we there- 
fore have to also consider the emission at shorter wavelengths, 
in the mid-infrared range. 

We constructed the spectral energy distribution (SED) of 
the above described 1° radius region (which includes the en- 
tire CNC, but excludes the clouds around Gum 31) in the fol- 
lowing way: First, we extracted the total fluxes in this region 
from our Herschel maps and found values of 844 500 Jy in the 
70yum band, 847200 Jy in the 160//m band, 374700 Jy in the 
250 fim band, 156000 Jy in the 350 ^m band, and 56 000 Jy in 
the 500 //m band. For fluxes at shorter wavelengths, we use the 
M SX and IRAS fluxes in the 8 - lOO/vm range listed in Table 1 
in Smith & Brooks! d2007l) . which were derived for a very simi- 
lar are^^FlrrallyTwe also added the FIR and mm-band fluxes de- 
termined in the recent stud y bv|Salatino et al.1 (120121) . who also 
considered the same area as lSmifh & Brooks! (120071) . The result- 
ing global SED of the CNC is shown in Fig. [9] 

In order to estimate the total cloud mass from this SED, 
we repeated the radia tive transfer modeling as described in 
Preib isch et al.l (1201 lah . We emphasize that this radiative trans- 
fer modeling is not intended to be detailed and highly accurate, 
but just to see whether we can reproduce the general shape of the 
observed SED with reasonable assumptions about the mass and 
large-scale density distribution of the surrounding clouds. We 
do not intend to model the small-scale structure of the individual 
clouds; instead, we simply assume a central source of radiation 
surrounded by a spherical envelope of dust and gas. Although 

7 ISmith & Brooks! d2007h used a rectangular box aligned to the galac- 
tic coordinate system with a width of 1.91° (galactic longitude) and a 
height of 2.74° (galactic latitude), that includes almost all clouds associ- 
ated with the Carina Nebula, but excludes most of the emission around 
Gum 31. 




Wavelength [/xm] 

Fig. 9. Spectral energy distribution of the Carina Nebula com- 
plex. The thick squares show the fluxes derived from our 
Herschel maps for a 1 ° radius area. The crosses show the mid- 
and far-i nfrared fluxes determin ed from the IRAS and MSX 
maps by Smit h & Brooksl d2007l) for a very similar region. The 
cross with the upward triangle shows the lower limit to the to- 
tal 870 fim flux derived from our LABOCA map that covered the 
central 1 .2 5° x 1.25° area. The di amond symbols show the fluxes 
taken from lSalatino et ail (120121) . The thick solid line is the sum 
of the spectrum resulting from the radiative transfer model and a 

power- law spectrum of the form F v — 1100 Jy ( 3 4 / m ) (fol- 
lowing lSalatino et al.ll2012l) . 



this is obviously a strong simplification, it provides the advan- 
tage that the temperature distribution of the gas is computed in a 
more physically meaningful way than adding up a few discrete 
graybody components. 

Since a full des cription of the radiativ e transfer modeling 
was already given in lPreibisch et al. li EoTTah . we summarize here 
just the main parameters and results: 

The inner and outer edges of the spatial grid were set at 
r = 10 18 cm (0.3 pc) and 41 pc. Follow ing the census of mas- 
sive stars in the CNC by 'Smith (2006), we assumed a total 
stellar luminosity of 2.4 x 10 7 L o and a typical temperature of 
T e ff = 44700 K. The model shown in Fig. [9] assumes that the 
density increases with distance from the center according to 
p(r) oc r L35 for radii up to 20.4 pc and then stays constant up 
to the outer radius at 41 pc; this should approximately match the 
conditions in the CNC, where most of the cloud material has al- 
ready been dispersed from the central region and is now located 
at the periphery of the complex. 

The model shown in Fig.|9]contains a total (dust + gas) mass 
of 889 140 M Q . The spectrum of this simple model matches the 
observed FIR fluxes for A > 100 pm quite well. The fluxes at 
shorter wavelengths are not so well reproduced, but the model 
fluxes generally are within a factor of ~ 2 - 3 of the observed 
values. We note that a good fit of this part on the SED is not 
expected from our simple model. The apparent discrepancy is 
probably the consequence of the fact that the true temperature 
distribution is much more complex than our simple model. Since 
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our model SED predicts systematically too high fluxes in the 
15 - 60 /jm range, the derived mass should be considered to be 
an upper bound to the total mass. 

We also note that the LABOCA 870 fan flux of 729 Jy is 
clearly just a lower limit, because the LABOCA map covers just 
the central ~ 50% of the analyzed region in our Herschel map. 
This also explains why the mass derived from this modeling is 
larger than our previous estimate 280 000 M Q ) based on the 
smaller area of the LABOCA map. 

It is interesting to compare our result to previous mass esti- 
mates of the CNC based on fit s of the SED with a few blackbody 
curves. The SED modeling of Smith & Brooks (2007), who used 
a combination of three discrete blackbody curves with temper- 
atures of 35 K, 80 K, and 220 K suggested a total (dust + 
gas) mass of =s 962 000 M G , slightly higher but still well con- 
sistent with our resu lt of < 890 000 M Q . The SED model of 
Salat ino et al.l d2012l) . which was based on a single temperature 
component with T = 34.5 K, yielded a total dust + gas mass 
of » 950 000 M Q . Despite the uncertainties (e.g., due to the dust 
opacity), the results of these different approaches to measure the 
total mass agree remarkably well. 

8. Conclusions and summary 

Our Herschel maps show the FIR morphology of the clouds in 
the CNC at unprecedented sensitivity and angular resolution. 
They reveal the very complex and filamentary structure of the 
clouds, which seems to be dominated by the radiation and wind 
feedback from the numerous massive stars. 

The total (gas + dust) mass in the CNC traced by Herschel is 
~ 650 000 M Q , and our radiative transfer modeling suggests that, 
including also warmer gas that is not well traced by Herschel, the 
total mass may be as high as < 890 000 M Q . Most of this mass 
resides in clouds with rather low densities. Nevertheless, there 
are about 105 000 M & in clouds above the extinction threshold 
of Ay > 3 and about 22 600M o above Ay > 7. The lower of 
these two extinction values (Ay > 3) is generally considered 
to constitute the threshold above which the gas is confined in 
gravitationally bound entities, making it available for future star 
formation. The Ay > 7 threshold is though to distinguish the 
gas that is dense enoug h to be directly in volved in active star 
formation. We note that Lad a et al.l d2010b have demonstrated a 
good correlation between the gas mass above A y > 7 and the star 
formation rate, that seems to be valid over a very wide range of 
cloud masses. This relation predicts a star formation rate of the 
order 10 3 M yr _1 for the CNC, which is well consistent with 
the integrated mass and age distribution of the young stars tha t 
have formed recently in these clouds (see Preibis ch et al.l201 Id) . 

Another important aspect is the ratio between molecular and 
atomic gas in the complex. Our Herschel maps trace the FIR 
dust emission and thus our total (gas + dust) mass estimate of 
~ 650000 - 890 000 M includes both components, i.e. molec- 
ular and atomic gas. The amount of the molecular gas can be in- 
ferred from molecular line observations. lYonekura et all fe 005) 
presented a wide-field 12 CO (J=l-0) survey of the CNC and de- 
termined the total gas masses in different sub-regions. Our 1 de- 
gree region we used for our CNC mass estimates includes their 
regions No 1, 2, 3, 7, as well as some parts of their region 4 (see 
their Table 1 and Fig. 3). Adding up their numbers for the total 
gas mass based on the 12 CO signal gives Mnco ~ 200 000 M . 
This comparison suggests that about 25% of all the gas in the re- 
gion is in molecular form, while 75% is in atomic form. We note 
that the quoted percentages suffer from considerable uncertain- 
ties: on the one hand, the CO data may seriously under-estimate 



the total molecular mass, and on the other hand, uncertainties of 
the assumed dust opacity and calibration problems may affect 
the estimate of the total gas mass. Nevertheless, the suggested 
high mass fraction of atomic gas may be interpreted as the conse- 
quence of the high level of massive star feedback that irradiates, 
heats, and disperses the clouds since several million years. The 
feedback has apparently transformed much, presumably most, of 
the original molecular cloud mass into atomic gas. 

The Herschel maps provide much more detailed informa- 
tion than reported in this paper. The Herschel fluxes have been 
used in our recent study of the infrared detected jet-d riving pro- 
tostars in the Carina Nebula (Ohlend orf et al.l 12012b . In forth- 
coming publications we will report results about the detailed 
small-scale structure of the clouds, the numerous point-like FIR 
sources (most of which are expected to be embedded protostars), 
as well as the relation between the CNC and the Gum 3 1 nebula. 
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